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ABSTRACT 

Simple theoretical arguments indicate that cooled interstellar gas in bright elliptical galaxies 
forms into a young stellar population having a bottom-heavy, but optically luminous IMF ex- 
tending to ~ 2 Mq. When the colors and spectral features of this young population are combined 
with those of the underlying old stellar population, the apparent ages are significantly reduced, 
similar to the relatively young apparent ages observed in many ellipticals. Galactic mergers are 
not required to resupply young stars. The sensitivity of continuous star formation to Lb and 
Lx/Lb is likely to account for the observed spread in apparent ages among elliptical galaxies. 
Local star formation is accompanied by enhanced stellar H/3 equivalent widths, stronger optical 
emission lines, more thermal X-ray emission and lower apparent temperatures in the hot gas. 
The young stars should cause M/ L to vary with galactic radius, perturbing the fundamental 
plane of the old stars alone. 



Subject headings: galaxies: elliptical and lenticular 
galaxies: interstellar medium - x-rays: galaxies 

1. INTRODUCTION 

Traditionally, elliptical galaxies have been re- 
garded as ancient stellar systems in which evo- 
lutionary processes have been exhausted or com- 
pletely arrested. However, recent observational 
and theoretical developments have led to a re- 
assessment of both the ages of the stars in ellipti- 
cals and the ages since the stellar system merged 
into its current configuration. Our objective here 
is to illustrate the possibly important contribution 
to the apparent global stellar age in ellipticals due 
to a population of young, intermediate mass stars 
formed from cooled interstellar (cooling flow) gas. 
We also review the related structural history of 
ellipticals. The high degree of structural regular- 
ity among ellipticals has been used to argue that 
they were among the first galaxies that formed. 
Of most interest are the remarkable thinness of 
the fundamental plane (Djorgovski & Davies 1987; 
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Dressier et al. 1987; Renzini & Ciotti 1993), the 
tightness of the color magnitude relation (Bower, 
Lucey & Ellis 1992) and the tight correlation be- 
tween central Mg2 line strength and central ve- 
locity dispersion a (Bender, Burstein, & Faber 
1993; Ziegler & Bender 1997), indicating an in- 
timate connection between parameters that char- 
acterize chemical and dynamical evolution. The 
small scatter about these relations is difficult to 
reproduce if there is a wide range of ages among 
ellipticals. Moreover, distant elliptical galaxies at 
redshifts z > 1 exhibit these same strong correla- 
tions (Aragon-Salamanca et al. 1993; Ellis et al. 
1997; Stanford, Eisenhardt & Dickinson 1998; van 
Dokkum et al. 1998b; Barger et al. 1998; Broad- 
hurst & Bouwens 1999), but are brighter overall 
(van Dokkum & Franx 1996), suggesting an early 
formation epoch z >2 for an fim — 0.3, Ha = 0.7 
universe (van Dokkum et al. 1998a). 



Set against this conventional picture are a num- 
ber of recent studies that suggest that most or 
many elhpticals have formed continuously over 
time by mergers that introduce morphological dis- 
turbances (Schweizer & Seitzer 1992) and new 
generations of young stars (Kauffmann & Char- 
lot 1998). Some observations of stellar photo- 
metric indices (Prugniel, Golev, & Maubon 1999) 
and deficiencies of E galaxies at redshifts z < 1 
(Kauffmann, Chariot & White 1996) support sig- 
nificant ongoing evolution, particularly for field el- 
lipticals, but other studies at these redshifts indi- 
cate a nearly constant E/SO space density (Im et 
al. 1996: Franceschini et al. 1998) and old forma- 
tion ages (Bernardi et al. 1998). 

The age of stars in ellipticals has also been 
questioned, but progress has been confounded by 
the age-metallicity conspiracy: youthful, metal- 
rich and old, metal-poor populations are nearly 
indistinguishable (Worthey 1997). However, Lick 
observers have demonstrated that Balmer absorp- 
tion lines can partially break this age-metallicity 
degeneracy (Gonzalez 1993; Faber et al. 1995; 
Trager 1997). By comparing the H/3 equivalent 
width (index) and a Mg -t- Fe photometric in- 
dex with expectations from Worthey's evolution- 
ary models, Gonzalez found a wide range of appar- 
ent ages in a sample of 40 (mostly field) ellipticals 
for which the mean age is only 8 ± 3 Gyrs. In ad- 
dition, stars in the inner regions of these galaxies 
{r ^ ''e/8) are about 3 Gyrs younger (and more 
metal-rich) than those at larger radii. These stel- 
lar ages are inconsistent with 13-15 Gyrs, formerly 
thought to be more likely. The age spread is larger 
for ellipticals of lower luminosity. Most of the H/3 
equivalent width is contributed by F and G dwarf 
stars near main sequence turnoff. 

2. CONTINUOUS STAR FORMATION 

In a series of recent papers we have presented 
detailed models of the evolution of elliptical galax- 
ies with an emphasis on the gas dynamics of the 
hot interstellar gas (e.g. Brighenti & Mathews 
1999a; 1999b). Galactic stars are assumed to form 
at t^s ~ 1 Gyr and the de Vaucoulcurs structure of 
the large ellipticals is constructed at <* = 2 Gyrs. 
Our calculations successfully reproduce currently 
observed interstellar (cooling flow) density, tem- 
perature and (iron) abundance profiles in massive 



ellipticals with a minimum of adjustable parame- 
ters. 

The global rate that interstellar gas cools can 
be estimated by dividing the observed (bolomet- 
ric) X-ray luminosity L^ by the enthalpy per gram 
of the gas, M = (2^.mp/hkT)L^ « 2.5 Mq yr"! 
where T — 1.3 x 10^ K is typical of large ellipti- 
cals {nip ~ proton mass; ^ — 0.62 = molecular 
weight). The total mass of gas that cools over 
cosmic time, several 10^" Mq, is only about 4 - 
5 percent of the total baryonic mass currently in 
stars. 

Two of the most perplexing and long standing 
problems concerning galactic cooling flows are (1) 
to determine where cooling to low temperatures 
actually occurs in the galaxies and (2) to deter- 
mine the final physical disposition of the cooled 
gas. The dropout or cooling of interstellar gas 
must occur over a substantial volume of the in- 
ner galaxy, but the radial mass profile of cooled 
gas cannot be predicted from first principles since 
it depends critically on entropy fluctuations ac- 
quired during a variety of complex processes (stel- 
lar mass loss, supernovae explosions, magnetic 
field variations, etc.). To accommodate this un- 
certainty, we have considered a variety of cool- 
ing dropout models in which the hot gas is as- 
sumed to cool at a rate {dp/dt)do = —q{r)p/tdo, 
where tdo — 5mpkT/2fipA is the local (constant 
pressure) cooling time and q(r) is an adjustable 
dropout function (Brighenti & Mathews 1999b). 
We compare computed interstellar properties with 
currently (t„ = 13 Gyrs) observed interstellar gas 
in the luminous eUiptical NGC 4472. Most of 
these models are unacceptable because the radial 
distributions of X-ray surface brightness "Exir), 
gas density n{r) or temperature T{r) disagree 
with profiles observed in this galaxy. Among the 
models considered, the simple constant q{r) — 1 
model gave the best results, although the agree- 
ment with observed Sj;(r) was still not perfect; 
we consider this model again here to estimate the 
mass dropout in NGC 4472. 

Regarding the second perplexing and long 
standing problem, it has long been speculated that 
the end product of the cooled gas are low mass, 
non-luminous stars (Fabian, Nulsen & Canizares 
1982; Thomas 1986; Cowie & Binney 1988; Ved- 
der, Trester & Canizares 1988; Sarazin & Ashe 
1989; Ferland, Fabian & Johnstone 1994). We 



have recently reconsidered the star formation pro- 
cess in ehiptical galaxy cooling flows and have 
concluded that the mass of stars in the dropout 
stellar population probably extends to ^ 2 M©, 
i.e. the dropout population is optically luminous 
(Mathews & Brighenti 1999a). As cold gas col- 
lects at a cooling site, it becomes gravitationally 
unstable at this limiting mass, setting a firm upper 
mass limit m„ on the IMF for stars forming in the 
central regions of massive ellipticals. The upper 
mass limit on the (bottom-heavy) IMF increases 
only modestly with galactic radius {^ 4 Mq at 
r ~ Tg) and is almost independent of time dur- 
ing the evolution of the cooling flow for redshifts 
^< 1. 

Additional support for the formation of opti- 
cally luminous stars in galactic cooling flows is 
provided by the thinness of the fundamental plane. 
For agreement with observed X-ray surface bright- 
ness distributions in ellipticals, most of the mass 
of cooled gas is concentrated well within r,, where 
the relatively small dropout mass can contribute 
substantially to the central mass and mass to 
light ratio determined from stellar velocities. If 
the dropout stellar population is assumed to be 
non-luminous, Mathews & Brighenti (1999b) have 
shown that the variation of dark dropout mass 
among ellipticals causes large, undesirable shifts in 
the fundamental plane that are incompatible with 
its observed thinness. However, these perturba- 
tions on the fundamental plane may be lessened 
or removed if the dropout stars arc luminous. 

3. Lb, B-V AND H(3 

Star formation in ellipticals is efficient in the 
sense that the total mass of HII gas and cold neu- 
tral or molecular gas at any time are both much 
less than the total mass of gas that has cooled. 
Therefore, the star formation rate {'^sfr) in NGC 
4472 is equal to the instantaneous rate that hot in- 
terstellar gas cools by radiative losses. The total 
accumulated mass that has cooled in NGC 4472 
since i, = 2 Gyr, Mdo{t), and ^sfa = dMdo/dt 
are shown in Figure la; these are based on the 
q = 1 model that best fits the X-ray observa- 
tions of NGC 4472 (Brighenti & Mathews 1999b). 
The mass dropout Mdo{r,tn) at t„ = 13 Gyrs for 
this model occurs mostly in r < r^ and the total 
dropout mass is Mdo{tn) = 4.7 x 10^" Mq. This 



is much less than the total current stellar mass in 
NGC 4472, M^t = 7.26 x 10" Mq, determined 
with M/Lb = 9.2. 

For the purpose of illustration, we assume that 
the old galactic stars can be approximated as a sin- 
gle burst stellar population having a Salpeter IMF 
from ?Ti£ = 0.1 to niy: = 125 Mq. By contrast, 
the younger dropout stellar population with vari- 
able SFR, "^ sFR{t), is assumed to have a Salpeter 
IMF from m^ = 0.1 to to„ = 2.5 Mq. Both of 
these IMFs arc available in single burst format 
in the 1999 Bruzual-Charlot library for isochronc 
synthesis spectral evolution (Chariot & Bruzual 
1991: Bruzual & Chariot 1993), assuming solar 
abundance. The total B-band luminosity of the 
dropout stellar population is 
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where iB.do{t) is the single burst B-band lumi- 
nosity per unit solar mass. The luminosity of 
the background old single burst population is 
LB,oid{t) = fmM^t^B,oid{t) where /„ is a coef- 
ficient of order unity that must be adjusted for 
agreement with the observed B-band luminosity 
of NGC 4472 (see below). 

The evolution of B-band luminosities and B-V 
colors are illustrated in Figure lb and Ic for each 
stellar population and for their combined radia- 
tion. For reference, the redshift z = 1 is shown 
at time t = 6.19 Gyrs (assuming Ho — 65 km 
s~^, flm = 0.3 and Qa = 0.7). We have chosen 
fm = 1.35 so that the total B-band luminosity 
LB,tot — Ls^oid + Ls^do = 7.89 x 10 Lb,q, ap- 
propriate for NGC 4472 at distance d = 17 Mpc. 
Although the dropout population currently con- 
tributes about 15 percent of the total B-band light 
{LB,do{tn) = 1.2 X 10^" Lb,q), its fractional con- 
tribution to the galactic light is quite constant for 
redshifts z < 1. The combined do-|-old population 
is only slightly bluer (by 6{B — V) ^ 0.03) than 
the old population and this difference is essentially 
constant for z < 1. Formally, the B — V oi the 
combined population indicates an age '^ 8.5 Gyrs 
that is less than that of the old population alone, 
12 Gyrs, but metallicity variations could produce 
a similar color variation. 

To estimate the equivalent width of H/3 we as- 
sume that the line width is similar for both popu- 



lations: 
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Here fg{R) is the ratio of hght from the old to 
dropout population within projected radius R nor- 
malized to the total ratio of old to dropout light; 
generally we assume /g = 1, corresponding to 
viewing the total light from both populations. 

As an illustration, we use single burst ewf3{t) 
from the 1999 Bruzual-Charlot (BC99) tables ap- 
propriate to the IMF of each population. As 
shown in Figure 2 the dropout population (with 
/g = 1) reduces the apparent age of the old 
population by '--^ 5 Gyrs, i.e. EWf3,tot{tn) ~ 
EWf3^oid{tn ~ 5 Gyrs), in agreement with observa- 
tions that correlate bluer colors with stronger H/3 
(Forbes & Ponman 1999). For smaller jg — 1/4, 
corresponding to viewing NGC 4472 within Tg , the 
apparent age is reduced by ^ 8.5 Gyrs. Actual 
observations of galactic cores view a fraction of 
both populations, i.e. the apparent H/3 age is 
aperture-dependent. For ra^ > 2, ewf3^do{t) and 
EWp^totitn) are insensitive to m„, i.e. ewp^doit) = 
s-w fj ^oid{t) can be assumed. We also used the pop- 
ulation code available at the Worthey website to 
determine both EWi3,tot{tn) and {B — V)tot{tn) for 
a Salpetcr dropout IMF from mi = 0.2 to ruu = 10 
M© . For this model the II/3 and {B — V) ages are 
8.5 and 9.5 Gyrs with fg = 1. These age uncer- 
tainties are consistent with {^ 35 %) errors inher- 
ent to population synthesis procedures (Chariot, 
Worthey & Bressan 1996; Worthey 1996). 

Clearly, however, the contribution of cooling 
dropout stars to the spectra of elliptical galaxies 
can explain the relatively young ages inferred from 
H/3 observed in some ellipticals even if the under- 
lying stars are very old. 

4. FURTHER DISCUSSION 

All massive ellipticals contain cooling interstel- 
lar gas so no comparisons with gas-free galaxies 
can be made. However, since the cooling flow 



mass dropout is centrally concentrated within Tg, 
the equivalent width of H/3, EWp, should increase 
toward galactic centers, in agreement with the ob- 
servations of Gonzalez (1993). Dropout star for- 
mation should be accompanied by an ensemble of 
additional observations at small galactic radii: H/3 
in emission from cooling clouds, enhanced X-ray 
surface brightness due to dense, locally cooling 
regions, and lower apparent X-ray temperatures 
which (with hydrostatic equilibrium) indicate in- 
terior masses less than the known stellar mass. 
Such mass discrepancies within ^ Q-lr,, are appar- 
ent in X-ray observations of bright Virgo ellipti- 
cals (Brighenti & Mathews 1997a). Some fraction 
of the total optical light {^ 15 percent in Figure 
1) in E galaxies comes from the dropout popu- 
lation with {M/LB)do « 4 < {M/LB)oid- There- 
fore, {M/Lb) should vary with galactic radius and 
produce a shift away from the fundamental plane 
defined by the old stars alone. 

The observed apparent H/3 age spread among 
ellipticals is large: 2 to 12 Gyrs in the Gonzalez 
(1993) sample, 5 to 12 Gyrs in the Fornax cluster 
(Kuntschner & Davies 1998) and 8 to 12 Gyrs in 
the Coma cluster (Jorgensen 1999). Such varia- 
tions can be expected if the cooling dropout pro- 
file q{r) differs among otherwise similar galaxies, 
represented here with the factor fg. But some of 
the age variation may arise from comparing ellip- 
ticals of greatly different Lb- Interstellar X-ray 
emission from ellipticals with Lb ^ 3 x 10^" is 
masked by stellar X-rays, but cooling flows still 
exist in these faster rotating, low Lb ellipticals. 
We have shown (Brighenti & Mathews 1997b) that 
large cold gas disks may form from cooling flow 
gas in low Lb ellipticals similar to the HI disks 
observed by Oosterloo, Morganti & Sadler (1999). 
Occasional star formation with normal IMFs and 
maximal EW/j (de Jong & Davies 1997) is there- 
fore expected in low Lb ellipticals. In luminous 
ellipticals (Lb ^ 3 x 10^°) m„ and EWi3 wiU gen- 
erally be lower due to isolated star formation in 
the high pressure ISM. But L^/Lb varies enor- 
mously among bright ellipticals of similar Lb and 
H/3 ages should reflect this same variation since 
the fraction of mass in the dropout population is 
proportional to Lx/Lb- As expected, ellipticals 
in Gonzalez' sample with Lb > 3 x 10^" and low 
L^/Lb - NGC 4649, NGC 7619 and NGC 7626 - 
also appear to be very old {^ 13 Gyrs). However, 



this interpretation is not entirely straightforward; 
elhpticals with large Lx/ Lb also have larger inter- 
stellar pressure which may result in to„ < 2 and 
lower H/3 indices. Although H/3 ages are insen- 
sitive to the IMF of the dropout population for 
m„ > 2.5 Mq, as m„ approaches ^ Q.SMq, EWfj 
decreases; this is just the range in m„ anticipated 
from our model of star formation in luminous el- 
lipticals (Mathews & Brighenti 1999a). The H/3 
index may not vary monotonically with L^/ Lb- 

Although the star formation process described 
here obviates the need for continued galactic merg- 
ing to account for the H/3 equivalent widths ob- 
served, we do not claim that recent merging in el- 
lipticals is non-existent or unimportant. We do 
note, however, that rather few images of ellip- 
ticals at small redshift indicate ongoing mergers 
with gas-rich, star- forming galaxies. But if such 
regular mergers do generate young stellar popula- 
tions in ellipticals, it may be possible to detect az- 
imuthal asymmetries in the stellar H/3 that reflect 
the orbital plane(s) of the newly-introduced stars. 
Young stars formed from cooling flow dropout 
are expected to be symmetrically disposed in the 
galactic potential, but their orbits may be more 
radial with narrower lines than those of the old 
stellar population. 

We are indebted to Stephane Chariot for pro- 
viding H/3 indices for the two populations con- 
sidered and to Stuart Norton for helpful advice. 
Studies of the evolution of hot gas in ellipti- 
cal galaxies at UC Santa Cruz are supported by 
NASA grant NAG 5-3060 and NSF grant AST- 
9802994 for which we are very grateful. FB is 
supported in part by Grant MURST-Cofin 98. 
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Fig. 1. — (a) Solid line: Cumulative mass of cooled 
(dropped-out) gas in NGC 4472 in units of 10^° 
Mq] Dashed line: Star formation rate {Mq yr^^) 
from cooled interstellar gas in NGC 4472. (b) 
B-band luminosity in 10^*^ Lb,q for the dropout 
stellar population [short dashed line), the old stel- 
lar population (long dashed line) and both popu- 
lations combined {solid line), (c) {B — V) colors 
of the old stellar population [long dashed line) and 
both populations combined (solid line). The dot- 
ted line is drawn parallel to the horizontal axis. 
Redshifts z — and 1 are indicated for a cos- 
mology having Ho = 65 km s^^, ilm — 0.3 and 
Qa = 0.7. 

Fig. 2. — H/3 equivalent widths in Angstroms for 
the dropout stellar population (short dashed line), 
the old stellar population (long dashed line) and 
both populations combined (solid line). The dot- 
ted line is drawn parallel to the horizontal axis. 
Redshifts are labeled as in Figure 1. 
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